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Abstract. MIPAS observations of temperature, water va- 
por, and ozone in October 2009 as derived with the scien- 
tific level-2 processor run by Karlsruhe Institute of Tech- 
nology (KIT), Institute for Meteorology and Climate Re- 
search (IMK) and CSIC, Instituto de Astrofisica de An- 
dalucia (IAA) and retrieved from version 4.67 level-lb data 
have been compared to co-located field campaign observa- 
tions obtained during the MOHAVE-2009 campaign at the 
Table Mountain Facility near Pasadena, California in Octo- 
ber 2009. The MIPAS measurements were validated regard- 
ing any potential biases of the profiles, and with respect to 
their precision estimates. The MOHAVE-2009 measurement 
campaign provided measurements of atmospheric profiles of 
temperature, water vapor/relative humidity, and ozone from 
the ground to the mesosphere by a suite of instruments in- 
cluding radiosondes, ozonesondes, frost point hygrometers, 
lidars, microwave radiometers and Fourier transform infra- 
red (FTIR) spectrometers. For MIPAS temperatures (version 
V4O_T_204), no significant bias was detected in the mid- 
dle stratosphere; between 22 km and the tropopause MIPAS 


temperatures were found to be biased low by up to 2 K, while 
below the tropopause, they were found to be too high by the 
same amount. These findings confirm earlier comparisons of 
MIPAS temperatures to ECMWF data which revealed simi- 
lar differences. Above 12 km up to 45 km, MIPAS water va- 
por (version V4O_H2O_203) is well within 10% of the data 
of all correlative instruments. The well-known dry bias of 
MIPAS water vapor above 50 km due to neglect of non-LTE 
effects in the current retrievals has been confirmed. Some 
instruments indicate that MIPAS water vapor might be bi- 
ased high by 20 to 40 % around 10 km (or 5 km below the 
tropopause), but a consistent picture from all comparisons 
could not be derived. MIPAS ozone (version V4Q 03 202) 
has a high bias of up to +0.9ppmv around 37 km which 
is due to a non-identified continuum like radiance contri- 
bution. No further significant biases have been detected. 
Cross-comparison to co-located observations of other satel- 
lite instruments (Aura/MLS, ACE-FTS, AIRS) is provided 
as well. 
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1 Introduction 

Altitude-resolved satellite measurements of atmospheric 
temperature, water vapor content and ozone mixing ratios are 
essential to obtain a global picture of the state of the atmo- 
sphere in the light of global change. One instrument provid- 
ing such data is the Michelson Interferometer for Passive At- 
mospheric Sounding (MIPAS) (Fischer et ah, 2008) onboard 
the Envisat research satellite. MIPAS is a mid-infrared limb 
emission Fourier transform spectrometer designed for global 
vertical profile measurement of temperature and many atmo- 
spheric trace constituents relevant to atmospheric chemistry 
and climate change. The measurement range of MIPAS ex- 
tends from the upper troposphere to the lower thermosphere. 
MIPAS temperature measurements are a target result in their 
own right, because global altitude-resolved temperature in- 
formation particularly in the upper stratosphere and above 
is limited. Beyond this, precise knowledge of temperatures 
is an essential precondition to trace gas retrievals, because 
the thermal emission of trace molecules depends strongly on 
temperature, and any temperature retrieval error will propa- 
gate onto the retrieved concentration profiles. Water vapor 
and ozone, also part of the MIPAS data product, are essential 
climate variables, contribute to the greenhouse effect of the 
atmosphere, are involved in atmospheric chemistry, and are 
tracers of atmospheric transport. 

There exist multiple processors for analysis of MIPAS 
spectra; this paper focuses on temperature, water vapor and 
ozone profiles retrieved with the data processor operated by 
the Institute of Meteorology and Climate Research (IMK) at 
the Karlsruhe Institute of Technology (KIT) in cooperation 
with the Instituto de Astrofisica de Andalucia (von Clarmann 
et ah, 2003b), which supports analysis of a greater variety 
of atmospheric species than the operational ESA processor 
(Ridolfi et ah, 2000; Raspollini et ah, 2006) and more obser- 
vation modes with an extended altitude range. From sum- 
mer 2002 to spring 2004, MIPAS measured in its original 
measurement mode at a spectral resolution of 0.025 cm -1 
unapodized. Retrieval of temperature, water vapor and ozone 
profiles has been described by von Clarmann et ah (2003b), 
Milz et ah (2005), and Glatthor et ah (2005, 2006), respec- 
tively. The functionality of the retrieval processor and the un- 
derlying radiative transfer code KOPRA (Stiller et ah, 2002) 
were validated by von Clarmann et ah (2003a, c, 2002), re- 
spectively. The resulting MIPAS data product was validated 
by Wang et ah (2004, 2005) for temperature, Milz et ah 
(2009) for water vapor and Steck et ah (2007) for ozone. 
After a failure of the interferometer slide in 2004, mea- 
surements at the original high spectral resolution were no 
longer possible, and from January 2005 on measurements in 
the new so-called optimized-resolution nominal observation 
mode were recorded at 0.0625 cm -1 unapodized. The re- 
trieval scheme had to be adjusted to the new measurement 
mode (von Clarmann et ah, 2009). This paper reports the 
first validation of these new data products, which took place 
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within the framework of the MOHAVE-2009 campaign at 
Table Mountain (California) in October 2009 (Leblanc et ah, 
2011), where a multitude of in situ, lidar and remote sensing 
instruments provided coincident measurements. 

2 MIPAS data and retrieval 

MIPAS on Envisat provides in its optimized-resolution nom- 
inal observation mode about 1300 radiance profiles per day, 
each consisting of 27 radiance spectra covering the altitude 
range of 6 to 70 km, and the spectral range of 4. 15 to 14.6 pm. 
The sun-synchronous orbit of Envisat at approx. 800 km al- 
titude allows coverage of the globe from pole to pole, with 
a horizontal sampling of 410 km along 14.4 orbits per day. 
The vertical sampling is 1.5 km up to 21 km altitude, 2 km 
up to 3 1 km altitude, 3 km up to 46 km altitude and 4 km 
above. The instantaneous vertical field-of-view covers 3 km, 
i.e. oversampling is achieved in the troposphere and lower 
stratosphere. Due to its emission sounding capability, MI- 
PAS records spectra of the atmosphere during day and night. 

Retrieval of temperature and trace gases from the 
optimized-resolution nominal observation mode at IMK/IAA 
is described in von Clarmann et ah (2009). The retrieval is 
based on constrained inverse modelling of limb radiances. 
The IMK/IAA processor performs regularized retrievals on 
a finer altitude grid ( 1 km gridwidth in the troposphere up to 
the middle stratosphere). Thus, stable solutions can only be 
obtained by regularization. While other MIPAS processors 
(Burgess et ah, 2004, 2006; Hoffmann et ah, 2008) regu- 
larize by the maximum a posteriori (also known as optimal 
estimation) method (Rodgers, 2000), the IMK/IAA proces- 
sor uses a smoothing constraint, which operates by weighted 
minimization of the squared first order finite differences of 
adjacent profile values, using a Tikhonov (1963) formal- 
ism. The intent of this choice is to make the resulting pro- 
files less dependent on the a priori profiles. For each target, 
dedicated spectral ranges, so-called microwindows, are used 
which were selected such that the total error consisting of 
measurement noise and parameter errors from the forward 
modeling is optimized. 

For the retrieval targets analysed in this paper, i.e. temper- 
ature, water vapor and ozone, a detailed description of the 
specific retrieval approach, microwindows and the estimated 
precision, accuracy and vertical resolution for the current 
data versions (version V4O_T_204, version V4O-H2O2203, 
and version V40 03 202) is given in von Clarmann et al. 
(2009). A summary of the relevant numbers, i.e. vertical res- 
olution, measurement noise error, total precision (including 
measurement noise and all parameter errors of random na- 
ture), total accuracy (including total precision and all sys- 
tematic error sources), and horizontal resolution along the 
line-of-sight (in terms of full width at half maximum of the 
horizontal averaging kernel) is provided in Table 1. 
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Table 1 . Data and error chararacterization of temperature, water vapor, and ozone retrieved from MIPAS level- lb version 4.67 spectra 
(optimized-resolution nominal observation mode) at IMK/IAA. 


Retrieval 

target 

Temperature 
(version V4O_T_204) 

Water vapor 
(version V4O_H2O_203) 

Ozone 

(version V4O_O3_202) 

Vertical resolution 

3.4 km (10 km) to 
1.9 km (40 km) 

2.3 km (20 km) to 
6.9 km (50 km) 

2.4 km (20 km) to 
3.5 km (50 km) 

Measurement noise 

0.2 K (10 km) to 
0.8 K (50 km) 

0.13 ppmv ( 1 0 km) to 
0.84 ppmv (50 km) 

0.03 ppmv (10 km) to 
0.08 ppmv (50 km) 

Total precision 

0.5 K (10 km) to 
1.4K (50 km) 

0.20 ppmv ( 1 0 km) to 
0.92 ppmv (50 km) 

0.07 ppmv (15 km) to 
0.28 ppmv (40 km) 

Total accuracy 

0.5 K (10 km) to 
2.1 K (50 km) 

0.34 ppmv ( 1 0 km) to 
1 .06 ppmv (35 km) 

0.07 ppmv (10 km) to 
0.78 ppmv (30 km) 

Horizontal resolution 

128 km (10 km) to 
402 km (40 km) 

206 km (10 km) to 
436 km (40 km) 

253 km (10 km) to 
405 km (40 km) 


Preliminary comparisons of the retrieved MIPAS temper- 
atures with ECMWF temperature fields indicated that there 
might be a systematic retrieval problem in the subtropics (25° 
to 40° N/S) below ~22 km: MIPAS temperatures seemed to 
be systematically higher by up to 2 K below the tropopause 
and lower by up to 2 K between the tropopause and ~22 km. 
Since any temperature retrieval error will propagate onto 
the retrieved concentration profiles, a careful validation of 
temperatures is most important. 

3 MOHAVE-2009 campaign 

The Measurements of Humidity in the Atmosphere and 
Validation Experiments (MOHAVE) 2009 campaign took 
place at the JPL Table Mountain Facility (TMF) at 34.4° N, 
1 17.7° W on 1 1-27 October 2009. MOHAVE-2009 was an 
extended version of the MOHAVE and MOHAVE-2 cam- 
paigns held at TMF in October 2006 and 2007. These cam- 
paigns, endorsed by the Network for the Detection of Atmo- 
spheric Composition Change (NDACC), allowed a thorough 
evaluation of the water vapor Raman lidar measurements up 
to the lower stratosphere by comparing to RS92 radiosonde 
and cryogenic frost point hygrometer profiles. 

Though lidar validation had again triggered the planning 
of the campaign, many other instruments and techniques 
joined the intercomparison efforts, leading to one of the most 
extensive atmospheric water vapor validation campaigns ever 
performed. The main goal of the campaign was to validate 
the water vapor measurements of four Raman lidars, two mi- 
crowave radiometers, two types of operational radiosondes, 
two types of frost point hygrometers, and an Infrared Fourier- 
Transform Spectrometer, as well as the column water mea- 
surements of a Ultra-Violet Fourier-Transform Spectrometer 
and two Global Positioning System (GPS) receivers. Mea- 
surements from five satellite instruments were included in 
the set of correlative data. Another goal of the campaign 
was to provide water vapor profiles from the ground to the 


mesopause without gaps. The third and last objective was to 
study water vapor variability in the UTLS in connection with 
the position of the subtropical jet near TMF. 

MOHAVE-2009 not only hosted all the instruments of the 
earlier campaigns in 2006 and 2007, but included three ad- 
ditional instruments and/or techniques, leading to the cor- 
relative measurement of temperature and water vapor from 
the ground to the mesopause, and ozone from the ground 
to the stratopause. To optimize the lidar range, the core 
of the campaign was centered near 19 October at the oc- 
currence of the new moon. Additional high priority nights 
(i.e. selected timing and increased density of the measure- 
ments and balloon launches) corresponded to the Aura Mi- 
crowave Limb Sounder (MLS), Aura Tropospheric Emission 
Spectrometer (TES), Aqua Atmospheric Infrared Sounder 
(AIRS), Atmospheric Chemistry Experiment (ACE), and 
MIPAS best coincidences near TMF. The campaign opera- 
tions were adjusted in real time following the most favor- 
able atmospheric conditions. High-resolution potential vor- 
ticity (PV) analysis and forecasts from the MIMOSA trans- 
port model (Hauchecorne et ah, 2002) supported the mea- 
surement planning. A more detailed description of the cam- 
paign operations and planning rationale is provided in the 
review paper by Leblanc et al. (2011). 

3.1 Operated instruments 

A detailed description of the measurement principles and in- 
struments operated during MOHAVE-2009 is provided in the 
review paper by Leblanc et al. (2011), and the dedicated 
articles in the present special issue on the MOHAVE-2009 
campaign (Hurst et ah, 2011b; Leblanc et ah, 2012a; White- 
man et ah, 2011). Here we give only a short introduction 
to the instruments which have been used within the valida- 
tion of MIPAS. A summary of the instruments’ main char- 
acteristics in terms of vertical resolution and measurement 
total uncertainty as retrieved from the data files is provided 
in Table 2. 
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Table 2. Summary of vertical resolutions und total uncertainties of coincident measurements. The total uncertainties are derived by averaging 
over all provided error profiles of a single type within 10 -km altitude bins, after having transformed the original provided measured quantity 
and its error into K or volume mixing ratio. The range of these mean errors over all altitude bins is given in the table. Uncertainties not 
provided in the data files but taken from references cited in the text are given in bold face. 


Insrument 

Measured quantity 

Vertical Resolution Uncertainty [K] or [%] 

TMW/TMF lidar 

H 2 0 [ppmv] 

150 m up to a 

5.2-30% 


O 3 strat [molm -3 

few kilometers 

11-28% 


O 3 trop [molm -3 ] 


5-13% 


T[ K] 


0.7-2 . 6 K 

STROZ lidar 

H 2 0 [g kg - 1 ] 

up to 1 .5 km 

9/8% 


O 3 [m -3 ] 


0.7-19% 


T[ K] 


0. 13-0.67 K 

ALVICE lidar 

H 2 0 [ppmv] 

up to 1 .2 km 

2.3-27% 


T[ K] 


0.18-0.33 K 

CFH 1 

H 2 0 [ppmv] 

meters 

7% 


O 3 [ppmv] 


5-10 % 


T[C] 


0.5 K 

NOAAFPH 1 

H 2 0 [ppmv] 

meters 

7% 


O 3 [ppmv] 


5-10 % 


T[ C] 


0.5 K 

RS92_GSFC 

T[C] 

meters 

0.5 K 

RS92JPL 

T [C] 

meters 

0.5 K 

WVMS 

H 2 0 [ppmv] 

12-15 km 

8% 

MIAWARA-C 

H 2 0 [ppmv] 

12-15 km 

12 - 21 % 

MkIV 

H 2 0 [ppmv] 

3-6 km 

5% 

Aura/MLS 

H 2 0 [ppmv] 

1 .5-3 km 

5.7-22% 


O 3 [ppmv] 


1.6-60% 


T[ K] 


0. 6-2.2 K 

ACE-FTS 

H 2 0 [ppmv] 

3 km 

5.1-25% 


O 3 [ppmv] 


2.5-27% 


T[ K] 


2 K 

AIRS 

H 2 0 [gkg -1 ] 2 

2-7 km 

18-21% 


T[ K] 


1. 0-2.6 K 


1 


T from radiosondes (RS92) and O 3 from ozonesondes flown with the frost point hygrometers. 2 


All profiles with errors larger than 30 % have been removed before comparison. 


3.1.1 Lidars 

The JPL water vapor Raman lidar at TMF (TMW) is a 
high-capability lidar system dedicated to the measurement 
of water vapor in the upper troposphere and lower strato- 
sphere (Leblanc et al., 2008, 2012a). The light emitted by a 
Nd:YAG laser at 355 nm is inelastically backscattered by at- 
mospheric nitrogen and water vapor, and collected at 387 nm 
and 407 nm respectively. After a few typical signal correc- 
tions (background correction (top end of the raw signals), 
saturation correction (non linearity at the bottom end of the 
signals), correction for molecular extinction along the laser 
beam path to the backscatter altitude and back), the ratio of 
the lidar signals collected in the water vapor and nitrogen 
channels is proportional to water vapor mixing ratio. The 
profiles are calibrated using external measurements, more 
specifically radiosondes during MOHAVE-2009. Systematic 
uncertainty is estimated to be 5-10% mainly depending on 


the calibration accuracy. Precision is mostly driven by ran- 
dom noise (photon counting), and typically ranges for several 
hours of integration from under 0.5 % in the mid-troposphere 
to 50 % in the UTLS. To mitigate this noise, the profiles are 
vertically smoothed and the resulting resolution ranges from 
150 m at the bottom to a few kilometers at 20 km. An in- 
depth description of the TMW can be found in Leblanc et al. 
(2012a). In addition, two other mobile lidar systems from 
the NASA-Goddard Space Flight Center (GSFC), referred to 
hereafter as “ALVICE” and “STROZ” lidars were employed 
during MOHAVE-2009 and used for MIPAS temperature, 
water vapor, and ozone validation. 

The ALVICE system (Atmospheric Laboratory for Vali- 
dation, Interagency Collaboration and Education) is a mo- 
bile facility that includes various atmospheric instruments in 
addition to the Raman lidar. The system provides, among 
other components, measurements of water vapor and rota- 
tional Raman temperature measurements, which were tested 
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for the first time during the MOHAVE-2009 campaign. The 
performance of the various components of the ALVICE sys- 
tem are discussed in Whiteman et al. (2011). The vertical 
resolution of the ALVICE system ranges up to 1.2 km in 
the upper parts of the profile. For the comparison to MIPAS 
measurements, we used temperature and water vapor mixing 
ratio measurements from ALVICE. For water vapor, the so- 
called best estimate profiles were used. This best estimate 
product merges the variably smoothed, 1 h sum and all night 
lidar profiles and includes a ground value of mixing ratio de- 
rived from ground-based in-situ sensors. The all-night lidar 
product includes a correction for signal-dependent bias be- 
lieved to be due to fluorescence of contaminants present in 
the lidar telescope. 

The Stratospheric Ozone (STROZ) lidar, operational since 
1989, was developed within GSFC Stratospheric Chemistry 
and Dynamics Branch to be an ozone, and temperature li- 
dar validation standard within NDACC (formerly NDSC) 
(McGee et al., 1991, 1995). Other measurement capabili- 
ties have been added over the years (aerosols in 1992, and 
water vapor in 2005). Currently the lidar transmits a pair 
of wavelengths, 308 nm from a XeCl laser and 355 nm from 
a high powered Nd-YAG laser. The STROZ lidar operated 
in three separate modes during MOFIAVE 2009. First, an 
ozone mode with a FOV of 2.3 mrad and transmitting at 
308 nm and 355 nm typically was used for two hours dur- 
ing which ozone, temperature, aerosol, and water vapor were 
retrieved. The second mode transmitted only 355 nm with 
the main telescope closed down to 1 .0 mrad, a mode during 
which aerosol, temperature, and water vapor was retrieved. 
The third mode consisted of transmitting only 355 nm with 
a FOV of 1.0 mrad, but with a filter, which blocked 355 nm 
from the first telescope while transmitting 387 and 407 nm 
radiation. The block was placed ahead of the collimation op- 
tics of the main telescope. No such filter was placed in the 
4" receiver linked to the second telescope. This mode re- 
turned only water vapor data. The filter was used because 
it was shown from a previous MOFIAVE campaign that flu- 
orescence excited by the 355 nm within the receiver chain, 
although small, can (and did in the STROZ lidar case) pro- 
duce a wet bias in the water vapor retrieval at high altitudes 
(low water vapor). The blocking filter greatly reduced but 
did not completely remove this interference from the STROZ 
data. The vertical resolution of the STROZ system ranges 
up to 1.5 km in the upper parts of the profile. A thorough 
description of the instrument and data will be presented in 
McGee et al. (2012). 

Two other lidars permanently deployed at TMF and op- 
erated by JPL acquired tropospheric ozone, stratospheric 
ozone, and middle atmospheric temperature profiles through- 
out the MOFIAVE-2009 campaign. The stratospheric ozone 
Differential Absorption Lidar (DIAL) is permanently de- 
ployed at TMF, and has acquired over 2000 routine profiles 
since 1989 (2h per night, 3 to 4 nights per week). It uses 
two XeCl Excimer lasers and one Nd:YAG laser to transmit 
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in the atmosphere at 308 nm (weakly absorbed by ozone) and 
355 nm (non-absorbed). The backscattered light is collected 
on a 0.91 m diameter Newtonian telescope and sent to an 
optical receiver where it is separated in 6 channels, includ- 
ing two pairs of Rayleigh channels (low and high intensity), 
and one pair of Nitrogen Raman channels (for the ozone and 
temperature retrieval in the lower part of the stratosphere of- 
ten affected by aerosols or thin clouds). Ozone is retrieved 
between the altitudes of 15 km and 50 km with a total uncer- 
tainty ranging from 5% (mid-stratosphere) to 15% (upper 
stratosphere), while temperature is retrieved between 15 km 
and 85 km with total uncertainty ranging from less than 1 K 
(lower stratosphere) to 5-10 K (75-85 km). A tropospheric 
ozone DIAL complements the stratospheric measurements 
since 1999. It uses a Nd:YAG laser with a fundamental at 
1064 mu quadrupled to 266 nm. The beam is splitted in two 
high-pressure cells filled with Flydrogen and Deuterium for a 
Raman-shift to 289 nm and 299 nm respectively. These trans- 
mitted UV wavelengths are in a more absorbing region than 
that used for the stratosphere in order to guarantee sufficient 
sensitivity in the ozone -poor troposphere. About 1500 tropo- 
spheric profiles (3-25 km with total uncertainty ranging from 
3 % to 15 %) have been acquired since 1999. 

3.1.2 Radiosondes 

Two types of meteorological radiosondes, designed for 
worldwide use on operational basis, were launched during 
MOHAVE-2009, namely the iMET-l-RSB and Vaisala RS92 
radiosondes. For validation of MIPAS, RS92 radiosonde data 
were used for temperature only, since the data on water vapor 
volume mixing ratio calculated from the relative humidities 
in the overlap region of MIPAS and RS92 measurements are 
not accurate enough for a meaningful validation. As shown 
by Flurst et al. (2011b) the iMet sonde temperature data are 
biased by 0.5 K versus the RS92 radiosondes, and the RS92- 
RS92 comparisons suggest that the total uncertainty in RS92 
temperature measurements is better than 0.5 K throughout 
the profile, which is consistent with an uncertainty analysis 
by Luers (1997) about an earlier version of this sensor. A to- 
tal of 58 RS92 radiosondes were launched during MOHAVE - 
2009. In 14 cases, two RS92 were mounted on the same 
balloon payload (“duals”). Data were received by two sep- 
arate ground systems, one (called RS92 JPL in the follow- 
ing) owned, launched, and processed by the JPL lidar group, 
and the other one owned and operated by the ALVICE group 
(called RS92_GFSC in the following). For the two systems, 
the processing software (digicora) version is slightly differ- 
ent, and the GSFC sondes include a GPS receiver, while the 
JPL ones do not. Although it is not mandatory to distinguish 
between the JPL and GSFC RS92s since the accuracies are 
equivalent, we have kept them separate in the following. Fur- 
ther measurements from radiosondes come from frost point 
hygrometer launches. 
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3.1.3 Frost point hygrometers 

The balloon-borne NOAA frost point hygrometer (FPH) was 
first flown over Boulder, CO, in 1980 (Oltmans et al., 2000) 
and, to date, has produced a 3 1 -yr record of stratospheric wa- 
ter vapor mixing ratios (Hurst et al., 2011a). Frost point hy- 
grometers measure the frost point of air passing through the 
hygrometer, from which the partial pressure of water vapor is 
directly calculated. The technique relies on the maintenance 
of a stable frost (ice) layer on a temperature-controlled mir- 
ror. Under equilibrium conditions, the ice surface temper- 
ature and water vapor content of the passing air are related 
through the Clausius-Clapeyron equation. Water vapor par- 
tial pressure is divided by the dry atmospheric pressure to 
yield the water vapor volume mixing ratio. 

Starting in 2003, the cryogenic frost point hygrometer 
(CFH) was developed in parallel to the FPH, with an em- 
phasis on reducing instrument size and weight and improv- 
ing frost layer stability (Vdmel et al., 2007). The CFH and 
NOAA FPH were developed independently, therefore there 
are subtle differences in the ways they operate. Neither in- 
strument requires water vapor calibration standards or a wa- 
ter vapor calibration scale; only the mirror thermistor must be 
calibrated with high accuracy and this is accomplished using 
NIST-traceable standards. 

Temperature and pressure measurements used to convert 
frost point hygrometer data into RH values and volume mix- 
ing ratios, respectively, are from the accompanying radioson- 
des on each balloon. Measurements of temperature and pres- 
sure from two different sondes are provided - one is an iMet- 
1-RSB sonde (also called iMet sonde in this paper), the other 
a RS92 radiosonde. We have used the temperatures from 
the RS92 sondes for comparisons to MIPAS. Because there 
is an altitude-dependent pressure bias between the two types 
of radiosondes, the water vapor mixing ratios calculated us- 
ing one set of pressure differs from the other. Above 20 km 
the pressure bias makes frost point hygrometer water vapor 
mixing ratios calculated using iMet sonde pressures 1-5 % 
higher than if calculated using RS92 sonde pressures (Hurst 
et al., 2011b). Below 20 km, however, pressure differences 
between iMet and RS92 sondes are smaller, substantially de- 
creasing the relative water vapor mixing ratio differences. 
We have used in our comparisons to MIPAS the water vapor 
mixing ratios calculated with RS92 sonde pressures. 

The ozone data for balloon flights are from ozonesondes 
flown with the frost point hygrometers. All the ozoneson- 
des were from the same manufacturer (EnSci) and were 
the same model (2Z). These ozonesondes are of the ECC 
(electrochemical concentration cell) type (Komhyr, 1969; 
Komhyr et al., 1995). There are no differences between the 
ozonesondes flown with NOAA FPH and CFH. 
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3.1.4 Microwave radiometers 

Two ground-based microwave radiometers participated to the 
campaign, namely the Water Vapor Millimeter-wave Spec- 
trometer (WVMS) permanently deployed by the US Naval 
Research Laboratory (NRL) at TMF (Nedoluha et al., 2011) 
and the portable Middle Atmosphere WAter vapor RAdiome- 
ter (MIAWARA-C) from the University of Bern, Switzerland 
(Straub et al., 2010). During a 5-month validation campaign 
the standard deviation of the MLS (version 2) -WVMS dif- 
ferences was shown to be 5 % from 26-70 km and the sys- 
tematic difference was within 8% throughout this altitude 
range. Both instruments use the pressure broadening of the 
water vapor rotational transition absorption line near 22 GHz 
for the retrieval of the altitude distribution of water vapor. 

The daily profiles during the MOHAVE-2009 campaign 
cover an altitude range from about 30 km (26 km for WVMS, 
and 33 km for MIAWARA-C) to 70 km with a vertical reso- 
lution of 12 to 15 km. The altitudes covered depend on the 
signal-to-noise ratios of the integrated spectra, which them- 
selves depend on the tropospheric conditions. Analysis of 
the MIAWARA-C forward and retrieval model provides an 
estimate of errors in the profiles which are typical for ground 
based 22-GHz water vapor radiometers. The total system- 
atic 2 -a error, taking uncertainties from the a priori temper- 
ature information, the calibration and the spectroscopy into 
account, is below 16% at all altitudes, while the random er- 
ror from measurement noise increases from 10 % at altitudes 
up to 50 km to 25 % between 50 and 70 km. 

3.1.5 FTIR ground-based spectrometer MkIV 

The MkIV FTIR spectrometer was designed and built at JPL 
in 1984 (Toon, 1991). Since then it has been operated on 
different platforms (ground-based, balloon-borne, and air- 
borne) in the framework of a large variety of different cam- 
paigns dedicated mainly to the investigation of stratospheric 
chemistry. The MkIV can measure high-resolution spec- 
tra (maximum optical path difference of up to 200 cm) and 
covers a very broad spectral range (650-5650 cm -1 ). For 
the MOHAVE-2009 campaign water vapor profiles were re- 
trieved following the method described in Schneider et al. 
(2010). The range of sensitivity for the MkIV instrument 
is limited from the ground to the upper troposphere which 
makes comparisons to MIPAS difficult due to a very small 
overlap range. 

3.2 Co-located satellite observations 
3.2.1 Aura/MLS 

Aura MLS was launched on 15 July 2004 into a near po- 
lar sun-synchronous orbit at 705 km altitude, with ascend- 
ing equatorial crossing time of 13:45 LT (Schoeberl et al., 
2006). It scans the Earth limb providing 240 scans per orbit. 
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spaced 165 km along the orbit track, and ~3500 vertical pro- 
files per day, with near pole-to-pole global latitudinal cover- 
age from 82° S to 82° N. MLS observes thermal microwave 
- far infrared emission from the Earth’s atmosphere in five 
spectral regions. Temperature is retrieved from the 118 GHz 
O 2 and 234 GHz 0 18 0 lines as described in Schwartz et al. 
(2008), while H 2 O is retrieved from measurements of the 
183 GHz H 2 O rotational line spectrum (Read et al., 2007; 
Lambert et al., 2007), and ozone is retrieved from the 236 and 
243 Ghz lines (Froidevaux et al., 2008). The MLS data pro- 
cessing algorithm is based on the optimal estimation method 
and uses a two-dimensional retrieval-approach to determine 
temperature, geopotential height and trace gas concentrations 
(Livesey et al., 2008). Most data products are retrieved on a 
fixed vertical pressure grid with 6 levels per decade change in 
pressure from the troposphere to the stratosphere. In the case 
of temperature and H 2 O (and ozone for data version 3.3), 
the vertical pressure grid is finer in the troposphere and the 
lower stratosphere, with 12 levels per decade change in pres- 
sure between 1000 and 22 hPa (0-25 km). For this study 
MLS version 3.3 (v3.3) data have been used, and geopoten- 
tial heights (GPH) provided within the data files have been 
used as altitude registration. This produces an altitude shift 
of 0 to 500 m over the altitude range of 0 to 55 km, which has 
been considered acceptable but should be kept in mind when 
analysing the comparisons. Read et al. (2007) and Lambert 
et al. (2007) have reported on vertical oscillations in v2.2 
H 2 O by up to 8% at 31.6 hPa which, however, have been 
eliminated in MLS version 3.3. 

MIPAS data recorded in the special Upper Tropo- 
sphere/Lower Stratosphere (UTLS-1) mode have already 
been compared to MLS v2.2 data by Chauhan et al. (2009). 
Temperature, water vapor and ozone were compared along 
latitudes, and on basis of global means. For temperature, ver- 
tical oscillations typically up to ±4 K were observed between 
MIPAS and MLS in the pressure/latitude range from 316.2 
to 100.0 hPa and 90° S to 90° N, respectively. MIPAS was 
colder by up to IK than MLS at the 21.5 hPa pressure level. 
In the middle and upper stratosphere MIPAS and MLS were 
biased within ±3 K. Differences up to ±4 K in the UTLS and 
in the stratosphere were in agreement with the differences ob- 
served for MLS in comparison to other satellite instruments 
(Schwartz et al., 2008). The global mean altitude-dependent 
bias between MIPAS and MLS temperatures varied between 
±2.5 K peak-to-peak and was within the combined system- 
atic errors over the complete pressure range 3 1 6-0. 1 hPa with 
small exceptions at 14.6 and 3.2 hPa. 

Regarding water vapor, MIPAS was wetter than MLS be- 
tween 316.2 and 177.8 hPa over the mid-latitudes and poles 
by about 50% (up to 100%), more prominently over the 
Southern Hemisphere. From 215.4 to 177.8 hPa over the 
sub-tropics and tropics, MIPAS H 2 O volume mixing ratios 
(VMRs) were drier by 10 % compared to MLS. In the lower 
stratosphere (146.7-56.2 hPa), the MIPAS/MLS water va- 
por comparison showed oscillations of ±1 ppmv (±10%) 
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over all latitudes. Oscillations up to 10% were also ob- 
served between 31.6-26.1 hPa for all latitudes. In the middle 
and upper stratosphere between 26. 1-0.2 hPa the agreement 
between MIPAS and MLS was within ±5 %. The oscilla- 
tions observed in the height range 31.6 to 26.1 hPa were due 
to a known artefact in MLS H 2 O v2.2 retrievals (Lambert 
et al., 2007) which has been removed for version 3.3. The 
altitude-dependent global mean water vapor bias was within 
±0.2 ppmv (±4%) between 100 to 0.1 hPa, with an excep- 
tion at 3 1 .6 hPa and 26. 1 hPa due to the problem in the MLS 
v2.2 data as mentioned above. From 316 hPa to 100 hPa the 
mean bias in water vapor varied from +12 % to —4 %. 

MIPAS ozone was up to 10% higher than MLS in the 
upper stratosphere (6.8-1.4hPa) over the southern mid- 
latitudes and subtropics. Over the tropics at 31.6 hPa a dif- 
ference (MIPAS minus MLS) of 0.5 ppmv or 20% was ob- 
served. Over the tropics and sub-tropics above the tropical 
tropopause layer (TTL) (100.0-68.1 hPa) relative differences 
as high as +90% to —80% were observed, while over the 
South pole in the height range 146.7 to 100.0 hPa differences 
were up to 10%. The high relative differences in the trop- 
ics and sub-tropics in the TTL were partly explained by low 
ozone concentrations and strong vertical gradients in this re- 
gion. The global altitude-dependent mean ozone bias was 
found to be 5% or less than 0.35 ppmv, with MIPAS being 
always higher than MLS, and peaking at about 5 hPa. 

3.2.2 ACE-FTS 

The Atmospheric Chemistry Experiment Fourier Transform 
Spectrometer (ACE-FTS) is the principal instrument on- 
board the Canadian SCISAT satellite (Bernath et al., 2005). 
SCISAT was launched on 12 August 2003 into a 74° in- 
clined orbit with an altitude of 650 km. ACE-FTS is a high- 
resolution FTS with the following specifications: spectra are 
recorded from 750 cm -1 to 4400 cm -1 (13.3 to 2.2 pm), at 
a resolution of 0.02 cm -1 (±25 cm maximum optical path 
difference). The instrument measures using solar occupa- 
tion and provides up to 30 measurements per day. It records 
one full spectrum in about 2 s with a signal-to-noise ratio be- 
tween 300:1 and 400:1 near the center of the wavenumber 
range. The delay between consecutive spectra gives a verti- 
cal spacing varying from 1.5 to 6 km depending on the an- 
gle between the orbit plane and the viewing direction with a 
maximum altitude resolution of approximately 3 km due to 
the field-of-view of the instrument (1.25 mrad). Details of 
ACE-FTS spectral inversion process are described in Boone 
et al. (2005). In a two-step process, temperature and pres- 
sure are retrieved from CO 2 transitions first and then these 
parameters are used retrieve the trace gas profiles. We have 
used the version 3 retrievals for validation of MIPAS data. 
This newest data version has reduced the occurrence of os- 
cillations in the temperature profiles and the microwindows 
for all trace gases have been updated. This dataset is in the 
process of being validated. 
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3.2.3 AIRS 

The Atmospheric InfraRed Sounder (AIRS) was launched 
into Earth orbit on 4 May 2002 on board the Aqua satel- 
lite, part of the NASA Earth Observing System (Chahine 
et ah, 2006). AIRS is a medium-resolution infrared grat- 
ing spectroradiometer. As a multi-aperture slit and pupil- 
imaging system, a diffraction grating disperses the incom- 
ing infrared radiation into 17 linear detector arrays compris- 
ing 2378 spectral samples. At long wavelengths, the spectral 
resolution is about 0.5 cm -1 decreasing to about 2 cm -1 at 
shorter wavelengths. The AIRS retrieval is based on iter- 
ative least squares physical inversion of clear column radi- 
ances following the approach of Chahine (1968, 1977). The 
retrieval of the AIRS water vapor profile uses a large set of 
channels associated with the strong 6 -pm water band, while 
temperature information is derived from the 15 and 4. 3 -pm 
CO 2 bands, and ozone is retrieved from the 9.6-pm ozone 
band (Susskind et ah, 2003, 2006). Water vapor amount is 
retrieved at twelve standard pressure levels between the sur- 
face and 100 hPa, though sensitivity is low for mixing ratios 
of about 10 ppmv or less (Gettelman et ah, 2004; Fetzer et ah, 
2008). AIRS water vapor retrievals have been validated ver- 
sus aircraft and balloon in situ measurements (Hagan et ah, 
2004; Gettelman et ah, 2004; Tobin et ah, 2006) and versus 
MLS (Fetzer et ah, 2008). 

4 Validation method 

The coincidence radius and time applied in this study were 
1000 km and 4h. If several coincident profiles to the 
same MIPAS profile or correlative measurement profile were 
found, we have used all coincident measurements, even if 
this introduces some interdependence in the data set. This 
was done in order not to reduce the size of the statistical en- 
semble. An overview of the numbers of coincidences is given 
in Table 3. Tests have shown that the conclusions from the 
comparison of all coincident measurements do not differ sig- 
nificantly from those where only unique MIPAS - reference 
pairs were used. 

Since most of the correlative measurements have a much 
different vertical sampling and resolution than the MI- 
PAS measurements, we have resampled the profiles x — 
{x\ ,...,x n ) T on a common altitude grid and degraded the 
better-resolved profile to the vertical resolution of the lower- 
resolved profile by application of the averaging kernel and a 
priori profile of the latter. Typically, profiles of lower vertical 
resolution are represented on a coarser altitude grid and vice 
versa. As a first step, both profiles are sampled on a common 
altitude grid. Resampling of a coarse profile x c on a fine grid 
can be written as 

-£ c f — Wx c , (1) 

where W is an interpolation matrix. The inverse operation, 
to map a high-resolved profile jcf on a less dense grid, is not 
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Table 3. Number of observations coincident with MIPAS for 
temperature, water vapor, and ozone. 


Instrument 

Temperature 

Water vapor 

Ozone 

TMW/TMF lidar 

22 

22 

27 

STROZ lidar 

31 

18 

27 

ALVICE lidar 

68 

74 

- 

CFH 1 

18 

18 

18 

NOAAFPH 1 

11 

11 

11 

RS92.GSFC 

44 

- 

- 

RS92JPL 

81 

- 

- 

WVMS 

- 

61 

- 

MIAWARA-C 

- 

116 

- 

MklV 

- 

454 

- 

Aura/MLS v3.3 

47 

43 

47 

ACE-FTS v3 

5 

5 

5 

AIRS 

2696 

1054 

- 


1 T from radiosondes (RS92) and O 3 from ozonesondes flown with the frost point 
hygrometers. 


unique but a reasonable recipe to achieve this is (Rodgers, 
2000, Chapter 10.3.1) 

*fc = Vx f , (2) 

where 

v=(w 7 ’wr 1 w 7 \ (3) 

which satisfies VW = I, I = unity. The application of the av- 
eraging kernel A c of the low-resolved profile x c to the better- 
resolved profile JCf under consideration of the a priori profile 
jt a of the low-resolved retrieval can either be performed on 
the coarse altitude grid 

ifc = A c Vx f + (I — A c )x a (4) 

or on the fine altitude grid 

* f = WA c Vxf+W(I — A c )x a . (5) 

We have chosen the intercomparison on the coarse grid, ac- 
cording to Eq. (4). For most intercomparisons in this paper, 
particular those of MIPAS versus in situ measurements or li- 
dar profiles, the correlative measurements were resampled on 
the MIPAS vertical grid and degraded to the MIPAS resolu- 
tion. Exceptions are profiles from MIAWARA-C, WVMS, 
MklV, and AIRS whose vertical resolution is worse than that 
of MIPAS. In these cases, the profiles were interpolated to a 
grid which was the combination of both original grids, con- 
taining all grid points from the one and the other original 
grid, and the MIPAS profiles were degraded with the aver- 
aging kernel of the correlative measurement where available 
instead. 

These transformations also have to be applied to the re- 
lated covariance matrices S. The transformation of the error 
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covariance matrix Sf of the better resolved measurement on 
the finer grid onto the coarser grid is 

S fc =A c VSfV r A c r . (6) 

For hybrid cases, e.g. when the coarser resolved profiles are 
represented on a finer grid than that on which the better re- 
solved data are represented, or if one grid is finer in one alti- 
tude region but coarser in another, the tools discussed above 
can easily be combined to suit the particular application. 

In case of MIPAS water vapor data there is another com- 
plication which is that instead of mixing ratios the logarithms 
of water vapor mixing ratios are retrieved; also the averaging 
kernels and covariance matrices refer to the logarithms of the 
water vapor mixing ratios. 

The application of MIPAS averaging kernels to a better re- 
solved profile on the basis of the coarse-grid averaging kernel 
Ai nc of the logarithm of the water vapor mixing ratio then is 

Xf c = exp(Ai nc Vln(.x:f) + (I - Ai nc )ln(* a )) . (7) 

Also the covariance matrix of the fine-grid correlative mea- 
surement has to be transformed into the log-space before 
the logarithmic averaging kernels of MIPAS can be applied. 
Equation (6) becomes 

S]nfc = A lnc VS lnf V r Af nc , (8) 

where Si n f is calculated from the original covariance ma- 
trix in the linear domain, Sf, by generalized Gaussian error 
propagation as 


Slnf = 

(ii- 

., 0 \ 

Sf 

■ 

*l;f 

0 > 


V 0 ’ • 

1 

x n;f ) 


v o, . 

1 

x n; f / 


The back-transformation of the error covariance matrix 
into the linear domain after application of the logarithmic 
averaging kernel is calculated as 




0 \ 



.. 0 \ 

Sf c = 



Slnfc 




l 0, . 

• • %n\fc ) 


^ 0, . 

• • %n\ fc / 


In the case when MIPAS water vapor is compared to pro- 
files from a measurement of lower resolution and coarser 
grid, Eq. (4) can be directly applied to MIPAS profiles in 
the volume mixing ratio (VMR) domain without any further 
complication. The transformation of the MIPAS logarithmic 
covariance matrix into the linear domain again applies the 
formalism of Eq. (10). This applies to the WVMS and the 
MIAWARA-C experiment. 

For AIRS, no averaging kernels are provided, and the data 
are provided on a vertical grid which is almost identical to the 
MIPAS grid. In this case, the data are compared as they are, 
without any transformation. Different altitude resolutions 
have to be kept in mind when the differences are explained. 


Mark IV water vapor retrievals are performed in the log- 
arithm domain, too. Their averaging kernels can be directly 
used to transform the logarithmic MIPAS covariance matrix: 

S] n fc = A, nc VSlnfV 7 A]^ c . (11) 

Table 4 presents a summary on the transformations applied 
to each MIPAS-validation data set pair. After these trans- 
formations of measurements and their error estimates to a 
common grid and after having degraded the better resolved 
profiles to the lower resolution of the other measurement, the 
comparison of data is performed. For evaluation of individ- 
ual pairs of correlative measurements x i;c and iq ; f c , we com- 
pare their differences to their combined accuracies (when- 
ever available; for some instruments only random error or 
measurement noise estimates are available which then are 
used instead) which are calculated as 


°7 ; diff 



+ cr 2 f 
i ;fc 


( 12 ) 


with <t 2 c being the variance of the coarser measurement and 

a l-fc being that of the degraded finer measurement on the 
coarser grid, both at altitude i. MIPAS error estimates in- 
clude measurement noise error, further random parameter er- 
rors and systematic errors. We have used the errors provided 
in the data files of the correlative instruments as total uncer- 
tainty without separating into various error sources. 

According to von Clarmann (2006), we first assess the 
bias between MIPAS and the correlative measurements, be- 
fore the precision validation is performed. The bias /;, 
is the mean difference between the MIPAS profiles and 
the coincident observations after convolution of the bet- 
ter resolved profile with the averaging kernel of the lower 
resolved measurement: 


^2n=l^ n J;c . i ; fc ) 

Ni ' 


bi = 


(13) 


Here the bias £>, is calculated independently for each alti- 
tude grid point i from the available N, coincident obser- 
vations. Nj can be different for different altitudes because 
the altitude coverage of a measurement system under assess- 
ment may vary from profile measurement to profile measure- 
ment. The standard error of the bias, which is also the bias- 
corrected root mean squares (rms) difference of the profiles, 
is calculated as: 


'A: bias : 


j (Xn.i:c %n,i\ fc bi )“ 

Ni(Ni- 1) 


(14) 


We consider the bias bj as clearly insignificant if the inter- 
val bi ±CT, ; bias includes zero. Additionally, we compare the 
bias to the systematic error of MIPAS (correctly it should 
be compared to the combined systematic error, i.e. square 
root of the sum of squared systematic errors, however, we 
take the error estimates of the correlative instruments as total 
(random) uncertainty and have not used any information on 
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Table 4. Transformations applied to each MIPAS-validation data set pair. 


Instrument 

Measured 

quantity 

Transformation 

applied 

TMW/TMF lidar 

T 

h 2 o 

0 3 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 
Eqs. (7), (8)— (10) on MIPAS grid, MIPAS Ai nc applied 
Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

STROZ lidar 

T 

h 2 o 

0 3 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 
Eqs. (7), (8)— (10) on MIPAS grid, MIPAS Ai nc applied 
Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

ALVICE lidar 

T 

h 2 o 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 
Eqs. (7), (8)— (10) on MIPAS grid, MIPAS Ai nc applied 

CFH 1 

T 

h 2 o 

0 3 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 
Eqs. (7), (8)— (10) on MIPAS grid, MIPAS Ai nc applied 
Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

NOAAFPH 1 

T 

h 2 o 

0 3 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 
Eqs. (7), (8)— (10) on MIPAS grid, MIPAS Ai nc applied 
Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

RS92_GSFC 

T 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

RS92JPL 

T 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

WVMS 

h 2 o 

Eqs. (4), (8)— (10) on union grid, WVMS A c applied to MIPAS 

MIAWARA-C 

h 2 o 

Eqs. (4), (8)— (10) on union grid, MIAWARA-C A c applied to MIPAS 

MklV 

h 2 o 

Eqs. (4), (1 1) on union grid, MklV A c applied to MIPAS ln[H 2 0] 

Aura/MLS 

T 

h 2 o 

0 3 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 
Eqs. (7), (8)— (10) on MIPAS grid, MIPAS Ai nc applied 
Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

ACE-FTS 

T 

h 2 o 

0 3 

Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 
Eqs. (7), (8)— (10) on MIPAS grid, MIPAS Ai nc applied 
Eqs. (4), (6) on MIPAS grid, MIPAS A c applied 

AIRS 

T 

h 2 o 

Interpolation to union grid, no averaging kernel applied 
Interpolation to union grid, no averaging kernel applied 


* T from radiosondes (RS92) and O3 from ozonesondes flown with the frost point hygrometers. 


the systematic error alone), in order to assess if the bias can 
be explained by known systematic uncertainties. 

The bias-corrected rms difference between coincident 
measurements <r/ ; diff is linked to the standard error of the 
bias by 

°7 ; diff — y/ N / tJ, ■ bias ■ (15) 

In the case of perfect co-incidences and valid random error 
estimates of both measurements, <T/ ; diff is expected to equal 
the combined single profile random error (see Eq. 12, but 
without consideration of systematic error terms) and thus is 
used for precision validation. 

For the standard approach, we have not separated the avail- 
able measurements into day and night profiles, although in 
the upper stratosphere and mesosphere, some effects of non- 
local thermodynamic equilibrium (non-LTE) triggered by il- 
lumination are known to be present in MIPAS profiles. These 
aspects are discussed in Sects. 5.1.5, 5.2.6, and 5.3.4, related 
to the assessment of systematic biases due to non-LTE ef- 
fects. Some instruments provide measurements from various 
integration intervals, for example 6 h versus 24 h measure- 


ments of the microwave instruments, or nightly means ver- 
sus 10-min measurements of TMF lidar. We have selected 
the nightly mean measurements from the lidars, the 24-h 
measurements from WVMS, and the 6-h measurements from 
MIAWARA-C for comparison, using the assigned measure- 
ment time in the file headers for determination of a potential 
coincidence. 

In all figures in the following, the differences provided are 
MIPAS profiles minus the correlative measurements, the one 
adjusted in vertical resolution by the averaging kernel of the 
other where appropriate, and brought to the same (coarser) 
vertical grid. For individual profiles, the blue curve repre- 
sents the MIPAS profile with the blue error bars representing 
the MIPAS error due to measurement noise, the green line 
and green error bars represent the original correlative mea- 
surement with its provided error, and the black line gives the 
correlative measurement transformed with the MIPAS aver- 
aging kernel. In case of comparison to WVMS, MIAWARA- 
C and MklV, the black line is the MIPAS profile transformed 
with the averaging kernel of those measurements. In case 
of AIRS, no degradation with the averaging kernels of AIRS 
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has been performed. In the right panel of all these figures, the 
difference of individual profiles is compared to the combined 
total errors of the two instruments according to Eq. (12). 

For averages over the coincident measurements, the blue 
and black line give the average of the MIPAS profiles and 
the averaging-kernel transformed correlative measurements, 
respectively, except for WVMS, MIAWARA-C, MkIV, and 
AIRS, in which cases MIPAS has been transformed. The 
bias is provided together with its standard error (shown as 
error bars) and the MIPAS systematic errors (dashed lines) 
in a second panel. In the third panel, the combined to- 
tal precision of individual measurement pairs according to 
Eq. (12), but without consideration of the MIPAS system- 
atic error (dashed lines), is compared to the bias-corrected 
rms differences (dotted lines). Again, in case of the correla- 
tive measurements we have not distinguished between vari- 
ous error sources, and the error as provided is used. For water 
vapor, profiles are presented on a logarithmic axis, and rel- 
ative differences are shown: these are the mean differences 
of the profiles given as percentage of the mean profile of the 
correlative measurement. 


5 Validation results 

5.1 Temperature 

5.1.1 Comparison to lidar temperatures 

Temperature measurements by lidars during the MOHAVE- 
2009 campaign are available from the instruments TMF li- 
dar, STROZ and ALVICE. While TMF lidar covers all al- 
titudes from the ground to the mesopause, STROZ mea- 
surements are available up to the lower mesosphere, and 
ALVICE measurements cover the troposphere and the 
lowermost stratosphere. 

Figure 1 shows the comparison of a MIPAS temperature 
profile measured on 1 8 October 2009 with the TMF lidar co- 
incidence of a nightly mean profile. The MIPAS - TMF lidar 
difference is mostly within or only slightly larger than the 
total error of MIPAS. Below the tropopause, MIPAS temper- 
atures are higher than TMF lidar temperatures by about 1 .5 to 
2 K, while just above the tropopause, they are lower by about 
the same amount, reproducing very well the well-known sig- 
natures found in differences relative to ECMWF data. Above 
the tropopause up to about 50 km, MIPAS is within 2 K of 
the TMF lidar temperature profile, while in the mesosphere, 
the differences are better than 3 K. The bias derived from all 
available 22 coincidences (see Fig. 2, top row) is less than 
1 K, except below 10 km, directly above the tropopause, 
near 42 km, and near 60 km. In the stratosphere, the bias 
is always negative, while in the troposphere, it is positive. 
The systematic errors of MIPAS cannot explain the bias (see 
Fig. 2, top row, middle panel). Further, the bias-corrected 
rms differences are far above the combined precision of both 
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Fig. 1. Left: single MIPAS temperature profile (blue) measured on 
18 October 2009 and a coincident profile (within 1000 km, 4 h) 
of nightly mean TMF lidar measurements (green). The solid black 
line gives the TMF lidar profile degraded with the MIPAS averag- 
ing kernel. The blue and green error bars are measurement noise 
errors of the MIPAS and the correlative measurement, respectively. 
At is the time difference in hours, Ad the spatial difference in kilo- 
meters, and Ala and Alo are the latitude and longitude differences, 
respectively, in degrees. Right: the black solid line provides the ab- 
solute difference between the coincident profiles, while the dashed 
line is the combined total error (including measurement noise, fur- 
ther random errors and systematic error components) of the two 
instruments, both on the coarser grid. For more details, see text. 

instruments (see Fig. 2, top row, right panel), which hints 
towards an underestimation of the random error of one or 
both instruments. Comparing the same set of TMF lidar ob- 
servations to ECMWF temperature profiles interpolated to 
the geolocations of the related MIPAS measurements (not 
shown) demonstrates that ECMWF is virtually bias-free to 
the TMF lidar measurements; this comparison hints towards 
a bias of —0.5 to — 1 K of MIPAS temperatures throughout 
the stratosphere. 

Comparison of MIPAS temperature profiles to STROZ li- 
dar measurements (see Fig. 2, middle row) indicates that no 
significant bias in MIPAS temperature data is present be- 
tween 18 and 50 km. MIPAS temperatures are higher by 
up to 2 K just below the tropopause, but lower than STROZ 
further down in the troposphere. Above the stratopause, the 
comparison indicates a strong low bias of MIPAS. The bias in 
the troposphere and mesosphere is much larger than the MI- 
PAS systematic errors; this means that the differences cannot 
be explained by known systematic uncertainties of MIPAS. 
The bias-corrected root mean squares differences are about 
twice as large as the combined precisions of the instruments. 

A comparison to ALVICE profiles is possible in the tro- 
posphere and lowermost stratosphere only. The already 
observed pattern of higher temperatures from MIPAS (1— 
2 K) below the tropopause and lower temperatures above the 
tropopause is also reproduced by the comparison to ALVICE 
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Fig. 2. Top row, left panel: average over all MIPAS temperature profiles (blue) for which coincidences with TMF lidar profiles have 
been found, together with the average of degraded TMF lidar profiles (black) which are coincident with the MIPAS profile. Top row, 
middle panel: averaged absolute differences together with their standard errors of the mean (error bars) and the combined systematic error 
components of the two measurements (dashed lines). Top row, right panel: bias-corrected root mean squares (rms) differences (dotted line) 
and the combined precisions of individual MIPAS and TMF lidar profiles (dashed line). Middle row: same as top row, but for STROZ lidar 
temperature profiles. Bottom row: same as top row, but for ALVICE lidar temperature profiles. 
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(see Fig. 2, bottom row). The MIPAS systematic errors and 
the combined precisions are much smaller than the bias and 
the bias-corrected rms differences, respectively. 

There is no clear explanation for the different behaviour of 
STROZ vs. TMF and ALVICE in the troposphere. STROZ 
data reveal a high bias of tropospheric temperatures com- 
pared to MIPAS, while TMF and ALVICE show a low bias. 
The ALVICE lidar is a different technique (Raman rota- 
tional) than STROZ and TMF, and is designed for tropo- 
spheric temperature recording. TMF and STROZ are de- 
signed for measurement of stratospheric and mesospheric 
temperature, and therefore are not so trustworthy well be- 
low the tropopause. Incomplete beam-telescope field-of- 
view overlap leads to temperature errors that can go both 
ways (i.e. too hot or too cold). Saturation effects may also 
impact the quality of the measurements. Finally, increasing 
aerosol optical thickness in the troposphere may be another 
source of error. 

5.1.2 Comparison to radiosonde temperatures 

The RS92 sondes flown together with the frost point hy- 
grometers and ozonesondes provide accurate measurements 
of temperature up to about 30 km. Comparison to single tem- 
perature profiles provided by the RS92 sondes flown together 
with the CFH frost point hygrometer (see Fig. 3, top row) 
show in general good agreement with MIPAS except some 
oscillations of the MIPAS temperature profiles with a period 
of ~5 km which are not present in the RS92 - frost point 
hygrometer profiles. The deviations between pairs of single 
profiles, however, are mostly larger than the estimated total 
error of MIPAS. The mean differences (see Fig. 4, top row) 
reproduce the already known high bias below the tropopause 
(~+l K within 5 km distance of the tropopause to +2.5 K 
below 10 km altitude) and low bias (1-2 K) just above the 
tropopause. The bias is considered significant below 15 and 
between 22 and 25 km and not explainable by known system- 
atic errors of MIPAS. The bias-corrected rms differences are 
much larger than the estimated combined precisions, hinting 
towards a severe underestimation of the random errors of one 
or both instruments, or deviations introduced by natural vari- 
abilities within the spacial and temporal coincidence ranges. 
Similar differences are found in the comparison to the tem- 
perature data provided from RS92 flown together with the 
NOAA FPH (see Fig. 3, bottom row, and 4, bottom row), 
however with tropospheric differences of +1 K to +3 K in- 
stead of +1 K to +2.5 K. A relative shift of the compared 
profiles in altitude by about 200 m would remove most of 
the differences, except the low bias of MIPAS directly above 
the tropopause. 

Besides on the frost point hygrometer balloon flights, 
RS92 sondes were also flown individually and provided 44 
(RS92_GSFC) and 81 (RS92JPL) temperature profile mea- 
surements up to 30 km altitude coincident to MIPAS obser- 
vations. A number of coincident single RS92 GSFC and 
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Fig. 3. Same as Fig. 1 but for RS92 radiosonde temperature profiles 
recorded during flights with the CFH (top) and NOAA FPH (bot- 
tom) frost point hygrometers. For a more detailed description see 
Fig. 1. 

RS92JPL profiles had physically unreasonable outliers in 
the upper part of the covered altitude range. These pro- 
files have been removed manually on basis of visual inspec- 
tion before calculating the mean differences. After removing 
the outlier profiles, the general picture of the RS92_GSFC / 
RS92 _JPL - MIPAS temperature comparison is the same as 
found in the other comparisons: the high bias in the tropo- 
sphere and low bias in the stratosphere of MIPAS tempera- 
tures is confirmed (see Fig. 6). The single differences oscil- 
late rather strongly with values of 0 to +3 K in the tropo- 
sphere and — 3 to 0 K in the stratosphere (see Fig. 5). Again, 
the bias-corrected rms differences are much larger than the 
estimated combined precisions, hinting at overly optimistic 
precision estimates or large natural variability. 

5.1.3 Comparison to temperatures from satellite 
instruments 

Within the period of the MOHAVE-2009 campaign coin- 
cidences with other satellite instruments providing temper- 
atures were found for Aura/MLS, ACE-FTS, and AIRS. 


www.atmos-meas-tech.net/5/289/2012/ 


Atmos. Meas. Tech., 5, 289-320, 2012 


302 


G. P. Stiller et al.: MIPAS validation by MOHAVE-2009 


Averages Bias 



Scatter 


-- comb, precision 
• rms , { 

0 2 4 6 8 

Diff. [K] 


Averages Bias Scatter 



Fig. 4. Same as Fig. 2, but for RS92 radiosonde temperature profiles recorded during flights with the CFFI (top row) and NOAA FPF1 (bottom 
row) frost point hygrometers. For a more detailed description see Fig. 2. 


Coincidences were only searched for within 1000 km around 
the Table Mountain Facility. The number of coincidences for 
Aura/MLS and ACE-FTS were rather sparse (see Table 3), 
partly due to differing local overpass times, partly due to 
differing observation geometries. Nevertheless we provide 
here the average differences in temperature for MIPAS versus 
these three instruments. 

The differences of MIPAS versus Aura/MLS (v3.3) mean 
temperatures provide a strong signature below the tropopause 
with MIPAS being warmer by up to 3 K in the troposphere; in 
contrast, MIPAS is colder than MLS in the stratosphere be- 
tween 25 and 42 km with the low bias peaking at 37 km and 
reaching —6 K (see Fig. 7, top row). The bias cannot be 
explained by the systematic error of MIPAS for most of the 
altitude range, and the estimated precison is about a factor 
of two to four lower than the bias-corrected rms, pointing 
towards underestimation of random errors. 

The MIPAS versus ACE-FTS mean temperature differ- 
ences oscillate within a band of ±3 K, with maximum devi- 
ations around 37 and 42 km (—3 K and +3 K, respectively), 
which, however, are not considered significant (the bias is 
not different from zero within its 1 -a uncertainty). Signif- 
icant biases are found between 19 and 23 km, and above 


50 km (see Fig. 7, middle row). The biases are larger than the 
estimated systematic errors of MIPAS. The bias-corrected 
rms differences are much larger than the estimated combined 
precisions of the two instruments up to the stratopause, but 
considerably smaller than the combined estimated precision 
above. The over-all structures of the two mean temperature 
profiles are consistent, although the stratopause seems to be 
lower by 2-3 km in the case of MIPAS which may partially 
explain the differences just below and above the stratopause. 
The AIRS mean temperature profile deviates from the MI- 
PAS mean profile in the troposphere (MIPAS being higher 
by up to 3 K), just above the tropopause (MIPAS being lower 
by 2 K), around 30 km (MIPAS being lower by 2 K), and in 
the stratopause (MIPAS being lower by 2 K). Over all the 
stratosphere, MIPAS tempera hues seem to be biased low ver- 
sus AIRS temperatures by about 1 K on average (see Fig. 7, 
bottom row), while in the troposphere, MIPAS is higher by 
about 2 K. The deviations are larger than the systematic er- 
rors of MIPAS, and the combined precisions are lower than 
the bias-corrected rms differences which hints towards an 
underestimation of the precision. 
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Fig. 5. Same as Fig. 1 but for RS92_GSFC (top) and RS92JPL 
(bottom) radiosonde temperature profiles. For a more detailed de- 
scription see Fig. 1. 

5.1.4 Non-LTE aspects 

The effect of non-LTE in the MIPAS temperature retrievals 
below 70 km over the Table Mountain Facility is not signifi- 
cant. This is because the population of the 15 pm states from 
which the temperature is derived, mainly that of the 0 1 1 0 vi- 
brational level, is in or very close to LTE at mid-latitudes 
below 85-90 km, even during daytime (Lopez-Puertas and 
Taylor, 2001). The temperature error caused by the non- 
inclusion of non-LTE is smaller than 0.2 K below 60 km and 
smaller than 0.5 K at 70 km. 

5.2 Water vapor 

Water vapor was the main validation target of the MOHAVE - 
2009 campaign. The goal of the MOHAVE-2009 cam- 
paign was to provide an accurate intercomparison of the in- 
struments widely applied to measure water vapor from the 
ground or from balloons. We took this validation oppor- 
tunity for comparison to the water vapor profiles derived 
from MIPAS. 
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5.2.1 Comparison to lidar water vapor measurements 

The TMW lidar provides water vapor measurements from 
ground up to about 22 km for nightly mean profiles, while 
MIPAS gives information from about 6 km (or cloud top al- 
titude) up to the lower mesosphere. An example for the 
comparison of a single MIPAS water vapor profile with a 
nightly mean TMW lidar observation on 17 October 2009 
is shown in Fig. 8. The differences between the single MI- 
PAS and the nightly mean TMW lidar water vapor profiles 
are within 10 % above 13 km, but reach +50 % below. Above 
13 km, the differences are smaller than the total estimated 
error of MIPAS. 

For the averages over all coincidences (see Fig. 9, top 
row), the differences between MIPAS and the nightly mean 
TMW lidar profiles are within 10%, except for the lower- 
most (below ~13 km) and uppermost (above ~24 km) alti- 
tude ranges, and especially the region around 10 km where 
the differences reach their maximum of +30 %. At 14 km and 
above, the bias between MIPAS and TMW lidar can fully 
be explained by the systematic errors of MIPAS which are 
driven by the uncertainties of spectroscopic data. The bias- 
corrected rms differences below 15 km are larger than the 
estimated combined precision which hints towards underes- 
timation of the random errors or a very high natural vari- 
ability; the latter often makes water vapor validation very 
difficult. In the stratosphere the bias-corrected rms differ- 
ences are even smaller than the combined precison which 
hints towards overestimation of the random uncertainty of 
the data and might be considered as another hint that natural 
variability may play a significant role below this altitude. 

The STROZ lidar covers water vapor up to about 
22 km; comparisons to MIPAS could be made up to 17 km. 
As reported in Leblanc et al. (2011) and Whiteman et al. 
(2011), STROZ lidar water vapor measurements at high al- 
titudes (low water vapor) are biased high due to undesired 
fluorescence, although a blocking filter was applied which 
greatly reduced, but did not completely remove the fluores- 
cence (see Sect. 3.1.1). Leblanc et al. (2011) reported that 
the wet bias of STROZ water vapor data started at 10 km and 
reached +20 % above 15 km. The MIPAS profiles are lower 
than STROZ by —30% and more below 10 km, and lower 
by —15 % to —20% at 12 km and above (see Fig. 9, middle 
row) the latter being consistent with the former findings. The 
bias-corrected root mean squares differences are by far larger 
than the estimated combined precisions. 

Similar to the STROZ water vapor profiles, the ALVICE 
water vapor profiles are restricted to altitudes below 
22 km; we compare here to the best estimate version of the 
profiles which has been corrected for a bias due to undesired 
fluorescence (see Sect. 3.1.1, Leblanc etal., 2011; Whiteman 
et al., 2011). The mean MIPAS water vapor profile above 
10 km agrees within 8% with the corrected, so-called best 
estimate mean ALVICE profile while below 10 km, MIPAS 
shows again a low bias. The bias-corrected rms differences 
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Fig. 6. Same as Fig. 2, but for RS92 J3SFC (top) and RS92 JPL (bottom) radiosonde temperature profiles. For a more detailed description 
see Fig. 2. 


agree roughly with the combined precision at the uppermost 
altitudes (above ~15 km) and are larger below (see Fig. 9, 
bottom row). 

5.2.2 Comparison to frost point hygrometer water 
vapor measurements 

The frost point hygrometers provide water vapor mixing ra- 
tio measurements of high precision and accuracy in the tro- 
posphere and lower stratosphere up to ~30 km. Deviations 
between CFFI and MIPAS individual water vapor profiles are 
sometimes very small and indicate that MIPAS reproduces 
the stmcture of the water vapor profiles, especially the posi- 
tion and deepness of the tropopause, very well (see Fig. 10, 
top row). The sharp structure with the sudden drop of water 
vapor VMR around 10 km, however, cannot be resolved by 
MIPAS, which is demonstrated by the CFH profile adjusted 
to the vertical resolution of MIPAS by applying its averaging 
kernel (black line), since the adjusted CFH profile is lower 
above and higher within the water vapor drop. But even com- 
pared to this degraded profile, MIPAS has a high bias below 


14 km which increases with decreasing altitudes. Compari- 
son to single NOAA FPH profiles confirm the comparison to 
CFH (see Fig. 10, bottom row). 

Figure 11, top row, shows the average over all MIPAS 
vs. CFH coincidences. Below appr. 13 km the differences 
exceed +20 %, while in the tropopause region and above the 
mean profiles agree within —5 and +8 %, and MIPAS repro- 
duces very well the profile shape. Above 13 km, the bias can 
fully be explained by the systematic errors of MIPAS, driven 
mainly by spectroscopic uncertainties. The bias-corrected 
rms differences and the estimated combined precisions are 
very close between 15 and 18 km and the estimated com- 
bined precision is even larger than the bias-corrected rms 
above, hinting towards a good error estimate of both instru- 
ments in this altitude range, and a too pessimistic precision 
estimate above 1 8 km. 

The comparison with the NOAA FPH reproduces in gen- 
eral the picture obtained from the CFH comparison, however, 
differences in the stratosphere are slightly higher and within 
— 12 and +5% (see Fig. 11, bottom row). The severe high 
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Fig. 7. Same as Fig. 2, but for Aura-MLS (top row), ACE-FTS (middle row), and AIRS (bottom row) temperature profiles. For a more 
detailed description see Fig. 2. 


bias of MIPAS below 10 km is confirmed. The bias-corrected 
rms differences are reasonably close to the estimated com- 
bined precisions around 15-18 km and higher above, again 
pointing towards a too pessimistic precision estimate in the 
stratosphere. 


5.2.3 Comparison to microwave radiometer water 
vapor measurements 

Microwave radiometers operated from the ground gener- 
ally provide a coarser altitude resolution in the stratosphere 
than MIPAS. For this reason, the MIPAS profiles have been 
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Fig. 8. Same as Fig. 1 but for TMW lidar water vapor profiles; the 
absolute water vapor VMRs are presented on a logarithmic scale, 
while the difference and the combined total errors are presented as 
percentage of the reference profile. For a more detailed description 
see Fig. 1. 

degraded with the averaging kernels of the microwave instru- 
ments to adjust vertical resolution. The WVMS instrument 
provides water vapor measurements from about 26 km to the 
mesosphere. In the stratosphere, the MIPAS average profile 
starts to deviate positively from the WVMS average profile 
around 30 km and develops an increasingly high bias which 
peaks around 37 km with a value of +10% (see Fig. 12, 
top row). A systematic low bias of WVMS in the strato- 
sphere on the order of — 10 % has been found in comparison 
to other instruments as well (Leblanc et ah, 2011) which is 
consistent with the comparison to MIPAS. Right above the 
stratopause, the difference between MIPAS and WVMS be- 
comes smaller, but a low bias of MIPAS in the order of — 8 % 
is found a few kilometers above in the mesosphere which can 
be explained by the neglect of non-LTE-effects in the MI- 
PAS retrievals (see Sect. 5.2.6). The deviations between the 
two average profiles are within the combined systematic er- 
rors everywhere. The bias-corrected rms differences and the 
combined precisions fit very well in the stratosphere, hinting 
towards a reliable error estimate of the instruments. 

The MIAWARA-C instrument provides water vapor pro- 
files from about 33 km to the mesosphere. MIPAS is lower 
than MIAWARA-C over all the comparison range, with a bias 
reaching —10%, except near 42 km where the bias is close 
to zero. Assuming an overall low bias of —10% between 
WVMS and MIAWARA-C, the high bias of MIPAS around 
40 km found in the comparison to WVMS is reproduced in 
relative terms by this comparison, and the low bias above 
50 km due to neglecting non-LTE is also confirmed (see 
Fig. 12, bottom row). Again, the bias can fully be explained 
by the combined systematic errors of the two instruments, but 
the combined precisions are about twice the bias-corrected 
rms differences which hints towards an overestimation of the 
random errors. 
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5.2.4 Comparison to water vapor profiles from the 
ground-based FTIR spectrometer MkIV 

Similar to the ground-based microwave instruments, the 
FTIR spectrometer MkIV operated from the ground provides 
a coarser altitude resolution than MIPAS, and MIPAS pro- 
files have been adjusted with the MkIV averaging kernels to 
allow meaningful comparison of the instruments. The sensi- 
tivity of MkIV reaches from the ground to the upper tropo- 
sphere (~15 km, as provided by the averaging kernel), and 
the overlap range where both MIPAS and MkIV are sensi- 
tive often is small. Contrary to most of the frost point hy- 
grometer and lidar instruments, the comparison of MIPAS 
to MkIV reveals a strong negative bias of MIPAS of up to 
—30% at altitudes below 15 km (see Fig. 13). Furthermore, 
the bias-corrected root mean squares differences are much 
larger than the combined precisions. This is inconsistent with 
the MOHAVE 2009 MkIV-RS92 comparison of Schneider 
et al. (2010), which showed good agreement between MkIV 
and the corrected Vaisala RS92 sondes. It should be kept in 
mind, however, that our formalism according to Eq. (4) and 
its variants disregards any a priori content of the better re- 
solved profile (in this case MIPAS), which might not be fully 
appropriate in the upper troposphere. 

5.2.5 Comparison to water vapor measurements from 
satellite instruments 

Comparison to other satellite instruments suffer from the 
few available coincidences found during the period of the 
MOHAVE-2009 campaign, and a satellite intercomparison 
would be better done globally. Nevertheless, we include here 
the coincidences found for October 2009 within 1000 km 
around Table Mountain Facility, to make intercomparison of 
all instruments possible. We have found 43 coincidences be- 
tween MIPAS and Aura/MLS which is mainly due to the 
different overpass times and a temporal coincidence crite- 
rion of 4 h (similar to the other comparisons). We have used 
Aura/MLS version 3.3 data for the comparison. Despite the 
sparse coincidences, the two instruments show good agree- 
ment, particularly in the stratosphere (see Fig. 14, top row) 
where the bias is between —6 and +10 % which can be fully 
explained by the combined systematic errors of the profiles. 
Below the hygropause the bias between MIPAS and MLS 
oscillates within ±20 %, while in the mesosphere, the well- 
known low bias of MIPAS is confirmed by the MLS ob- 
servations. The combined precision is larger then the bias- 
corrected rms above 13 km which hints towards an overly 
pessimistic random error estimate. 

Although coincidences between MIPAS and ACE-FTS are 
sparse, the comparison between the two instruments pro- 
vides a picture consistent with the other comparisons (see 
Fig. 14, middle row). The two instruments cover a similar 
altitude range from the upper troposphere to the lower meso- 
sphere. Average deviations between the two mean water 
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Fig. 9. Same as Fig. 2 but for water vapor profiles from TMW lidar (top row), STROZ lidar (middle row), and ALVICE lidar (bottom 
row); the absolute water vapor VMRs are presented on a logarithmic scale, while the bias and the various errors and bias-corrected rms 
differences are presented as percentage of the average reference profile. For a more detailed description see Fig. 2. 


vapor measurements are over all well below ±10% except 
below 15 km, where differences reach —20 to —25%, and 
above 60 km, caused by the non-inclusion of non-LTE in 
our MIPAS retrievals. It is to be noted that this is one of 
only few comparisons, besides MkIV and the lidars STROZ 


and ALVICE, where MIPAS is biased low in the troposphere. 
The estimated combined precisions again are higher than the 
bias-corrected rms differences and confirm the finding of an 
overly pessimistic precision estimate in the stratosphere. 
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Fig. 10. Same as Fig. 8 but for CFH (v3.20, top row) and NOAA 
FPFI (bottom row) frost point hygrometer profiles of water vapor. 
For a more detailed description see Fig. 8. 


AIRS retrieves water vapor between the surface and 
100 hPa, though sensitivity is low for mixing ratios of about 
10 ppmv or less (Gettelman et al., 2004; Fetzer et al., 2008). 
Since this is the lower part of the MIPAS observations, com- 
parison is somewhat difficult, despite the high number of co- 
incidences found. The bias between MIPAS and AIRS varies 
between ±10% and can be explained by the systematic er- 
rors of MIPAS, however, it is mostly not significant because 
of its large standard error (see Fig. 14, bottom row). The up- 
permost data point of the comparison is just at the edge of 
the sensitivity range of AIRS and should therefore be taken 
with some caution. Both the estimated combined precision 
and the bias-corrected nns are large, with the combined pre- 
cision being about twice the bias-corrected rms, which hints 
towards an overestimation of the random errors. 

5.2.6 Non-LTE aspects 

The effects of non-LTE on the water vapor retrieval at mid- 
latitudes are more important than on kinetic temperature. 
The daytime populations of the water vapor (010) vibrational 


level, from which water vapor abundance is retrieved, de- 
parts from LTE as low as 60 km. This is due to the strong 
coupling of that level with the 02 ( 1 ) level, which is pop- 
ulated after ozone photolysis. That produces water vapor 
(010) populations larger than in LTE. The strong water va- 
por fundamental band lines used in our retrievals are still at 
these altitudes under an optically thick regime. Instead of in- 
creasing the local water vapor abundance around 60 km, the 
global fit technique used in the water vapor retrieval com- 
pensates the smaller radiance simulated with the lower LTE 
populations at that tangent height by decreasing the water 
vapor abundance at altitudes above, reducing that way the 
absorption along the line of sight. Figure 15 shows the effect 
on the retrieved daytime water vapor at mid-latitudes if the 
non-LTE effects are included. The figure presents the mean 
difference of about 100 profiles retrieved from MIPAS nom- 
inal mode measurements with and without taking non-LTE 
into account in the retrieval. The MIPAS LTE retrievals un- 
derestimate water vapor by 20 % at 70 km, 8 % at 60 km and 
overestimate it by 5 % at 45-50 km. Opposite to the behavior 
above 60 km, the LTE retrieval at lower altitudes responds to 
the reduced absorption along the light of sight by increas- 
ing the abundance at that altitude. Although of systematic 
nature, the error due to neglect of non-LTE effects in the re- 
trieval cannot be easily corrected, since this correction would 
depend on illumination (non-LTE excitation), the actual ki- 
netic temperature profiles, and the actual amounts of other 
atmospheric constituents, and would thus not be constant for 
different atmospheric conditions. 

In summary, these non-LTE effects on the water vapor re- 
trievals can explain the differences found with other instru- 
ments in the lower mesosphere and also half of the difference 
found in the upper stratosphere. The rest of the difference in 
the upper stratosphere is in most cases anti-correlated with 
differences in the kinetic temperature, thus, the latter are their 
most likely source. 

5.3 Ozone 

5.3.1 Comparison to lidar ozone measurements 

The lidar instruments provide ozone measurements from the 
upper troposphere to about 45 km, while MIPAS measure- 
ments reach up to 70 km in the nominal observation mode. 
At TMF two lidars are operated optimized for measure- 
ment of tropospheric and stratospheric ozone, respectively. 
For more details, see Sect. 3.1.1 and Leblanc et al. (2011). 
The comparison of MIPAS mean ozone profiles to strato- 
spheric ozone measurements by the TMF lidars is shown 
in Fig. 16, top row. A peak in the MIPAS profile around 
37 km with positive deviations of 0.7 ppmv is obvious, while 
the other parts of the tropospheric and stratospheric profiles 
agree within 0.3 ppmv. The deviations between MIPAS and 
the TMF lidar are within the range of the combined system- 
atic errors. The combined precisions, however, are smaller 
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Fig. 11. Same as Fig. 9 but for CFH v3.20 (top row) and NOAA FPH (bottom row) frost point hygrometer profiles of water vapor. For a 
more detailed description see Fig. 9. 


than the bias-corrected rms differences above 25 km, but in 
agreement below. The tropospheric ozone measurements by 
the TMF lidars are shown in Fig. 16, middle row, together 
with the MIPAS profiles. The agreement is very good and 
does not exceed the combined systematic errors, while the 
bias-corrected rms differences are somewhat larger than the 
combined precisions. Both the stratospheric and the tropo- 
spheric TMF lidar ozone measurements hint towards an os- 
cillation in the MIPAS profiles with maximum values around 
22 and minimum values around 27 km with a peak-to-peak 
amplitude of 0.3 ppmv, which, however, does not exceed the 
estimated systematic errors of MIPAS. 

The STROZ lidar provides ozone measurements between 
approx. 15 and 45 km; comparison of the STROZ mean pro- 
files to MIPAS mean profiles is shown in Fig. 16, bottom row. 
The comparison is consistent with the findings from TMF li- 
dar: MIPAS ozone has a bulge with high ozone values around 
37 km which, however, is still within the estimated system- 
atic error. The positive bias of MIPAS reaches +0.9 ppmv at 
37 km and remains below ±0.3 ppmv for all other altitudes 
below 45 km. The oscillation in the difference profile with 
a maximum at 20 and a minimum at 27 km shows up in the 
comparison to STROZ as well, again with a peak-to-peak 


amplitude of about 0.3 ppmv. The bias-corrected rms dif- 
ferences are about twice as large as the estimated combined 
precisions above 25 km which hints towards underestimation 
of the random errors for one or both instruments. 

5.3.2 Comparison to ozonesonde measurements flown 
with the frost point hygrometers 

The ozone data provided with the frost point hygrometer 
measurements are from ozonesondes that were flown to- 
gether with the frost point hygrometers. There are no dif- 
ferences between the ozonesondes flown with NOAA FPFI 
and CFFI, so in principle the results from CFFI and NOAA 
FPFI balloon flights could be combined. We have kept them 
separate in order to follow the general scheme of compar- 
isons to all instruments. The ozonesondes flown with the 
frost point hygrometers provide ozone measurements below 
~30 km. As expected, the results from the comparisons 
to CFFI and NOAA FPFI are rather similar. MIPAS mean 
ozone profiles agree very well with the mean profiles of 
the ozonesondes, with an overall negligible bias, but an os- 
cillation in the difference profile with maximum at 20 and 
minimum at 24 km (see Fig. 17). The amplitudes of the 
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Fig. 12. Same as Fig. 9 but for WVMS (top row) and MIAWARA-C (bottom row) microwave radiometer water vapor profiles (6 h measure- 
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Fig. 13. Same as Fig. 9 but for MkIV FTIR spectrometer water vapor profiles. For a more detailed description see Fig. 9. 


oscillations are different: they remain below ±0.25 ppmv 
in case of CFH (Fig. 17, top row), and below +0.15 and 
—0.2 ppmv in case of NOAA FPFI (Fig. 17, bottom row) 
which might be due to differences in the actual spatial and 
temporal mis-matches between MIPAS measurements and 


the balloon flights. The estimated combined precisions are in 
both cases considerably smaller than the bias-corrected rms 
differences between 15 and 24 km, which is a hint towards 
underestimated random errors. 
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Fig. 14. Same as Fig. 9 but for water vapor profiles from Aura/MLS v3.3 (top row), ACE-FTS (middle row), and AIRS (bottom row). For a 
more detailed description see Fig. 9. 


5.3.3 Comparison to ozone profiles from satellite 
instruments 

From the satellite instruments, Aura/MLS, ACE-FTS, and 
AIRS provide ozone. The AIRS profiles, however, were 
found to be a scaled a priori profile only without providing 


information on the specific profile shape. For this reason, 
we compare to Aura/MLS and ACE-FTS only. The mean 
ozone profiles of Aura/MLS and MIPAS have rather dif- 
ferent shapes. They agree reasonably well, with deviations 
of less than ±0.3 ppmv, below 30 km, but have different 
shapes around the ozone VMR maximum; in particular, the 
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Fig. 15. Mean difference between MIPAS LTE and non-LTE water 
vapor daytime retrievals at mid-latitudes. About 100 profiles have 
been averaged. 

bulge in MIPAS profiles around 37 km shows up here, too 
(see Fig. 18, top row). Above the stratospheric ozone VMR 
maximum, the profiles are more or less parallel, but shifted 
in altitude by 3 km and more. Resulting deviations above 
30 km are in the order of —0.8 and +1.2 ppmv which is 
no longer explainable by the MIPAS systematic errors. The 
bias-corrected rms differences and the estimated combined 
precision are close below 25 km and above 35 km, while in 
between the bias-corrected rms differences are larger. 

The mean ozone profiles of MIPAS and ACE-FTS agree 
quite well in shape (see Fig. 18, bottom row); the differences 
reveal similar features as found in the comparison to lidar 
instruments: a positive bias of MIPAS peaking at 37 km and 
exceeding the systematic error estimate of MIPAS, a negative 
bias reaching —0.5 ppmv above, and some oscillating struc- 
tures, although less pronounced, around 20 and 24 km. ACE- 
FTS ozone is known to be biased high between 45 and 60 km, 
although reduced in version 3 from the version 2.2 ozone 
update, which explains the difference between MIPAS and 
ACE-FTS in this altitude range. The bias-corrected rms dif- 
ferences and the estimated combined errors agree well above 
39 km and below 22 km, while in between the precision 
error seems to be underestimated or the natural variability 
was large. 

5.3.4 Non-LTE aspects 

The ozone retrieval microwindows cover emissions from the 
i>2 and the V3 levels, in particular, from the fundamental and 
first hot bands. Their populations depart from LTE above 60- 
65 km during daytime, mainly due to the recombination of 
molecular and atomic oxygen, which produces vibrationally 
excited ozone. The ozone overestimation due to neglecting 
non-LTE is smaller than 1 % below 55 km and increases to 
20 % around 65-70 km (Gil-Lopez, 2006). Since MIPAS 
ozone profiles are lower than the correlative measurements in 


the mesosphere, neglect of non-LTE in the MIPAS retrievals 
cannot explain the differences. 


6 Conclusions 

MIPAS measurements of temperature, water vapor and ozone 
from the upper troposphere to the lower mesosphere retrieved 
from level- lb so-called optimized-resolution spectral data 
(version 4.67) with the IMK/IAA processor have been com- 
pared to balloon-borne and in-situ measurements performed 
during the MOHAVE -2009 campaign at Table Mountain Fa- 
cility, California in October 2009, and to co-located satellite 
instruments. The coincidence criteria were 1000 km in dis- 
tance and 4 h between the co-located measurements. All co- 
incidences between MIPAS profiles and correlative measure- 
ments were considered. We analysed both individual pairs of 
profiles and averages over all coincidences. We compared the 
mean difference profiles to their standard errors and the es- 
timated combined systematic error in order to detect signifi- 
cant biases not explained by known systematic error sources. 
Further we compared the bias-corrected rms differences to 
the estimated combined precisions in order to judge if the 
precision estimates were realistic. 

6.1 Synopsis from all instruments: temperature 

The temperature comparisons with most instruments pro- 
vided a consistent picture: in the stratosphere, no significant 
bias was detected. Both the altitude and the amplitude of the 
stratopause is well captured by MIPAS. Comparison with the 
TMF lidar and AIRS hint towards a low bias of MIPAS of 
1 K throughout the stratosphere, however, this bias is not 
significant over wide parts of the profile and was not con- 
firmed by other instruments. Comparison to MLS tempera- 
ture profiles do not agree with the consistent picture from the 
other instruments: MLS has a pronounced high bias versus 
MIPAS in the tropopause and is further biased high between 
25 and 42 km. Differences in temperatures may propagate to 
species retrievals: MIPAS and MLS water vapor differences 
seem to be anti-correlated to differences in their tempera- 
tures; temperature differences could also explain the larger 
differences in ozone found between MIPAS and MLS in the 
stratosphere. Around the tropopause, MIPAS has revealed a 
high bias on the order of 2 K below the tropopause, and a 
low bias of the same amount in the lowermost stratosphere. 
This behaviour has been suspected before from comparisons 
with ECMWF analysis data where it occured in the subtrop- 
ics only. Further down in the troposphere, most instruments 
indicate that the high bias of MIPAS remains between 1 and 
2 K, while STROZ and MLS indicate a low bias of about 2 K. 
In the mesosphere, no consistent picture could be gained, but 
there is a tendency of MIPAS temperatures being too low 
around 60 km. The mean difference profiles between MIPAS 
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Fig. 16. Same as Fig. 2, but for the stratospheric TMF lidar ozone profiles (nightly means) (top row), the tropospheric TMF lidar ozone 
profiles (nightly means) (middle row), and STROZ ozone profiles (nightly means) (bottom row). For a more detailed description see Fig. 2. 


and all other instruments, together with their standard errors, 
are compiled in Fig. 19, top panel. 

The detected bias profiles are in general larger than the 
estimated systematic error profiles, i.e. the bias cannot be ex- 
plained by known systematic uncertainties which are driven 
by spectroscopic uncertainties in case of MIPAS. The bias- 


corrected rms differences are typically between 2 and 3 K 
with a pronounced maximum around 17 km reaching val- 
ues up to 5 K. In particular the latter indicates that part of 
the bias-corrected root mean squares difference may come 
from high natural variability within the coincidence radii (in 
space and time). Leblanc et al. (2011) (their Fig. 8) and 
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Fig. 17. Same as Fig. 2, but for ozonesonde profiles recorded on flights with the CFH v3.20 (top row) and NOAA FPH (bottom row) frost 
point hygrometers. For a more detailed description see Fig. 2. 


Leblanc et al. (2012b) showed that during the MOHAVE - 
2009 campaign the TMF site was just at the edge of the 
subtropical jet stream, and a stratospheric intrusion was ob- 
served passing the TMF site on 20 October 2009, which 
both caused high natural variability in temperature, water 
vapor and ozone around the tropopause, even within a few 
hundred kilometers. 

The current analysis confirms earlier MIPAS temperature 
validation work of an older and different MIPAS tempera- 
ture data version by Wang et al. (2005), who also found a 
small overall bias in MIPAS temperature data. Similar to 
the current study, Wang et al. (2005) found also rather high 
bias-corrected root mean squares differences (2.5 to 3.5 K in 
their case); they assessed the contribution of natural variabil- 
ity to the overall bias-corrected root mean squares differences 
(rms) and found that more than 70% of the rms can be ex- 
plained by natural variability within the coincidence radii. In 
contrast to these earlier findings, however, is the high bias 
below the tropopause/low bias above the tropopause which 
showed up in the present assessment. 

The recent MIPAS temperatures (version V4O_T_204) 
were retrieved with a retrieval set-up which was different 
from the older one because it was adjusted to the lower 


spectral, but higher spatial resolution of MIPAS measure- 
ments since 2005. As a consequence of the comparisons to 
MOHAVE-2009 campaign data, we re-analysed the spectral 
ranges used within the temperature retrieval. This showed 
that one small spectral window used for the retrievals was 
contaminated by an ozone line, and was therefore sensitive 
to errors in the ozone climatology used. Test retrievals have 
demonstrated that the deviations around the tropopause al- 
most disappear if ozone is joint-fitted within the temperature 
retrieval so that it no longer depends on the used climatol- 
ogy. MIPAS temperature retrievals from version V40 T 205 
onwards will therefore be performed with the improved 
retrieval set-up including ozone as a joint-fit parameter. 

6.2 Synopsis from all instruments: water vapor 

Between 13 km and 55 km, the MIPAS water vapor mean 
profiles (see Fig. 19, middle panel) are within ±10 % of 
the profiles of the correlative measurements, except for the 
STROZ lidar, AIRS, and MkIV. The STROZ lidar instrument 
is known to have a high bias of up to 20 % above 12 km (see 
Sect. 3.1.1 and Leblanc et al., 2011) which explains the low 
bias of MIPAS versus STROZ. The discrepancy versus AIRS 
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Fig. 18. Same as Fig. 2, but for Aura/MLS (top row) and ACE-FTS (bottom row) ozone profiles. For a more detailed description see Fig. 2. 


is in the uppermost data point only which is around 100 hPa, 
and therefore at the edge of the sensitivity of the AIRS instru- 
ment; for this reason this discrepancy should be taken with 
some caution. The discrepancy versus MkIV is not in agree- 
ment with former MkIV validation efforts (Schneider et ah, 
2010). The discrepancy between TMW and MIPAS around 
25 km can be explained with the fact that TMW lidar data are 
valid up to 22 km only. The agreement with ACE-FTS and 
MLS is well within 10% all over the stratosphere. The mi- 
crowave instruments WVMS and MIAWARA-C both agree 
with MIPAS within ±10%. Three out of the four strato- 
spheric data sets point towards a tendency of MIPAS to be 
biased high up to 10% around 45 km, while above 55 km, 
MIPAS has a tendency to be too dry. Both biases are a well- 
explained consequence of the current MIPAS retrieval set up 
of water vapor which ignores non-LTE effects in radiative 
transfer in the mesosphere. The bulge at 45 km has been 
shown to be caused by error propagation from above and as 
such a compensation effect of the too low water vapor values 
retrieved in the mesosphere. In the troposphere below 13 km, 
no consistent picture could be achieved: while the frost point 
hygrometers and the TMW lidar point towards a high bias of 
MIPAS, STROZ and MkIV and, below 9 km, ALVICE and 
TMW indicate that MIPAS is biased low. 


The systematic errors of MIPAS, given by spectroscopic 
uncertainty and, above 40 km, non-LTE effects, can in most 
cases very well explain the biases to the other instruments. In 
the stratosphere, the bias-corrected rms is often lower than 
the estimated combined precisions by a factor of 2 which 
hints towards a far too pessimistic precision estimate; in the 
troposphere and around the tropopause, high natural vari- 
ability due to the vicinity to the subtropical jet stream may 
explain the larger rms compared to the precision estimates. 

An earlier MIPAS water vapor version (V30_H20_13) has 
been validated by Milz et al. (2009); they found no significant 
bias and a confirmation of the precision estimate which is on 
the order of 5-10 %. In particular, in their comparison to MI- 
AWARA measurements taken during a campaign in North- 
ern Finland, they found a similar bulge in the differences, 
but even more pronounced, around 45 km (see Milz et al., 
2009, their Fig. 16). Tropospheric water vapor data were not 
compared. The findings of the MOHAVE-2009 campaign are 
in good agreement with the previous findings by Milz et al. 
(2009), although the MIPAS observation mode was changed 
to lower spectral and higher spatial resolution in the mean- 
time. This is a good confirmation that the current retrieval 
setup is in accordance with previous data versions. 
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Fig. 19. Compilation of the mean differences between MIPAS and 
all instruments (MIPAS - correlative instrument) to which compar- 
isons have been performed. Top: all mean temperature differences 
and their standard errors; middle: all mean water vapor differences 
and their standard errors (in percent relative to the reference pro- 
file); bottom: all mean ozone differences and their standard errors. 
Vertical lines are meant as guide for the eyes only. 
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In the next data version of MIPAS IMK/IAA water vapor 
we will include non-LTE modelling in the radiative transfer 
calculations for the retrievals according to Garcia-Comas et 
al. (2012), which is expected to solve the problems at and 
above 45 km. 

6.3 Synopsis from all instruments: ozone 

A synthesis of all ozone comparisons is shown in Fig. 19, 
bottom panel. The comparisons to all relevant instruments 
provided the following picture: MIPAS ozone profiles have 
a pronounced high bias at the upper edge of the strato- 
spheric ozone VMR maximum around 37 km, with differ- 
ences reaching +0.9 ppmv in some cases (see Fig. 19, bottom 
panel). Between 50 and 60 km, the only instrument for com- 
parison is ACE-FTS. MIPAS is lower than ACE-FTS by up 
to —0.5 ppmv. However, a high bias of ACE-FTS ozone be- 
tween 45 and 60 km is a known feature, although reduced in 
version 3 from the version 2.2 ozone update. Below 30 km, 
the bias does not exceed ±0.3 ppmv except for the compari- 
son with MLS. In the lower stratosphere, an oscillation with 
a maximum around 20 km and a minimum around 27 km has 
been identified in several difference profiles, which, however, 
does not exceed the estimated systematic errors of MIPAS 
ozone, and does not show up consistently in all comparisons. 

A previous MIPAS ozone data version (V30_03_7) has 
been validated by Steck et al. (2007); in their comparisons, 
the bias between MIPAS ozone and other instruments was 
below ±0.3 ppmv, except for comparisons with HALOE and 
ground-based FTIR (see their Fig. 13). In particular, the 
mean comparison to lidars was better than 0.2 ppmv below 
40 km. Although not explicitly mentioned, the high bias 
found in the current data version around 37 km was present 
in some individual comparisons in version V30_03_7 as well 
(see their Fig. 5, top panel, or their Figs. 1 1 and 12). 

During analyses performed as a consequence of the find- 
ings from the present validation study, the ozone peak around 
37 km has been traced back to the handling of underly- 
ing continuum-like emissions in the spectral data. The 
continuum-like contribution is suspected to be caused by 
straylight from the Earth surface or the lowermost parts of 
the atmosphere, being scattered into the instrument’s optics 
by internal parts. In the current retrieval set-up, atmospheric 
continuum extinction and emission due to aerosols and other 
atmospheric constituents are accounted for by fitting an op- 
tical depth profile up to 32 km. A straylight-related radi- 
ance contribution at higher tangent altitudes can hence not 
be corrected. The straylight aspect is currently under fur- 
ther intense analysis; for MIPAS retrievals in the next future 
the continuum-like emission identified in the radiance spec- 
tra which leads to the ozone high bias around 37 km will be 
tackled as caused by an unidentified grey body which will be 
joint-fitted. Test retrievals have demonstrated that the high 
bias can be removed by extending the joint retrieval of an 
unidentified continuum extinction and emission up to 50 km. 
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